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The evolution of universe in Brans-Dicke (BD) theory is discussed in this paper. Considering a 
parameterized scenario for BD scalar field (f> — (f>oa°' which plays the role of gravitational constant G, 

o ■ 

' we apply the Markov Chain Monte Carlo method to investigate a global constraints on BD theory 

^ with a self-interacting potential according to the current observational data: the Union2 dataset of 

' type supernovae la (SNIa), the high-redshift Gamma- Ray Bursts (GRBs) data, the observational 

Hubble data (OHD), the cluster X-ray gas mass fraction, the baryon acoustic oscillation (BAO), 
and the cosmic microwave background (CMB) data. It is shown that an expanded universe from 
deceleration to acceleration is given in this theory, and the constraint results of dimensionless matter 
O . density f7om and parameter a are, ^om = 0.286l°;°39lo;o4? and a = 0.0046+o:oi7ilo o2og which is 

consistent with the result of current experiment exploration, | a |< 0.132124. In addition, we use the 
geometrical diagnostic method, jerk parameter j, to distinguish the BD theory and the cosmological 
constant model in Einstein's theory of general relativity. 
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I. Introduction 



lO . The observation of the supernovae of type la provides the evidence that the universe is undergoing accelerated 
expansion. Combining the observations from Cosmic Background Radiation and lager scale structure 4], one 
concludes that the universe at present is dominated by 70% exotic component. An interpretation to the accelerating 
universe can just be obtained by introducing an exotic fluid in the theory of general relativity. Considering the 
standard cosmological model this unknown component is called as dark energy with ovifning a character of negative 
pressure to push the universe to an accelerated expansion. And the popular dark energy model is the positive tiny 
cosmological constant (CC), though it suffers the so-called fine tuning and cosmic coincidence problems. However, in 
2(7 confidence level, it fits the observations very well Q] . If the cosmological constant is not a real constant but is time 
variable, the fine tuning and cosmic coincidence problems can be removed. In fact, this possibility was considered 



in the past years j6l-ll8l|. such as quintessence, Chaplygin gas, holographic, agegraphic dark energy, etc. Also, the 



accelerating universe is related to the modification of the gravity theory on large distances, such as f{R) modified 
gravity theory 19] and higher dimensional theory 20|, |21|, etc. Scalar-tensor modified gravity theories have recently 
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attracted much attention, in part because they emerge naturahy as the low-energy hmit of string theory as a result 
of the dilaton coupling with gravitons 22l-l27l|. In addition, it has been argued that in the early universe, gravity 



obey . s»,„ tensor tyE^ttafies „.her than general relativity (GR) M, and tlrey „e .l.oimportant tor 



cosmological inflation models 



23 



28 



30| . since they have effectively solved the problems of inflation [2^, where the 



end of inflation can be brought by nucleation without considering the fine-tuning cosmological parameters. 



II. Brans-Dicke theory 



Brans-Dicke (BD) theory 3ll433| of gravity is a simple but very important one among the scalar tensor theories, 

!„' i, 

which is apparently compatible with Mach's principle [34]. The generalized BD theory is an extension of the original 
BD theory by considering coupling parameter w as a function of the scalar field 35l438| , and an accelerating universe 
can be obtained when coupling parameter lo varies with time 39, 4^. In the framework of generalized BD theory, 
the action is described as (choosing the speed of light c — 1) 



S = d''x^\(j)R 



(1) 



where Lm is the matter Lagrangian, (j) is the BD scalar field which is non-minimally coupled to gravity, V{<j)) is the 
self-interacting potential for the BD scalar field, w(0) being a function of is generalized version of the dimensionless 
BD coupling parameter u, and <j) plays the role of the gravitational constant G which is related to the inverse of cj). 
The gravitational field equation derived from above action by varying the action with respect to the metric is 
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with 



□ 0: 



SttT 
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2(f> 



(2) 



(3) 



o-r2a;(0) 3-H2w((/))' ' d0 ^ 3 2w(<?!)) ^'''^ ' 

where G^i/ is the Einstein tensor, denotes the energy-momentum tensor of matter, and T — Considering 
Robertson- Walker universe 



-dt"" +a^{t)[- 



dr^ 



r^{d0^ 



,d^% 



(4) 



1 - kr^ 

where k denotes the spacial curvature with k ~ —1, 0, +1 corresponding to open, fiat and closed universe, respectively. 



The modified Friedmann equation and wave equation for the generalized BD scalar field are given as 
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43| 
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and 



(j) (j) 2(j) 



3 + 2a;(0) 3 2w((^) 
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[2V{4>) - 



d(f) ^ 3-f2w((/))^ 
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where pm and p denote the energy density and pressure for the matter, Pe// and Peff are effective energy density and 
pressure for the combination of matter and BD scalar field. Obviously, the conservation equation of effective fluid 
is satisfied, peff + 3H{peff + Speff) = 0. BD theory as a leading alternative to Einstein's theory of GR, it can be 
obtained by considering uj{4>) = uj =constant in above generalized form. Also, GR as a special case in BD theory, 
the theory can be approximately reduced to GR when take constant coupling parameter a; — >■ cxd and scalar function 
<f> ^constant. 



III. Cosmic constraints on Brans-Dicke theory 

It has been shown that the BD theory is effective for interpreting the c osmic acceleration by considering a non- 
minimal coupling between the gravitational term and the BD scalar field 
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-43|. In this work, we investigate the 
global cosmological constraints on BD theory. In order to solve above equations of motion the scale factor a = a^t^ is 
considered in keeping with the recent observations |44j|. Then the BD scalar field is solved as = as a function 
of time [4j|. Here /? and 7 are two constant parameters. Also, these two power-law solutions of the scale factor a and 



the BD scalar field 6 can be found in other references 



45l-l51l|. which are consistent with the evolution of expanded 



universe. Inspired by these papers, for calculation we consider a parameterized power-law form of Brans-Dicke scalar 
field, = 000." which can be easily given by above two solutions of the scale factor a and the scalar field (f). Then the 
Friedmann equation in BD theory with a flat geometry fc = is derived as, 
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(6 -|- 6a — a;a^)0o 
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(8) 



Defining ^ , ^ ^ = then one lias 



SttGq 0oF(0) 
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(9) 



with r^on 



In addition, for a flat universe according to Eqs. ([5]), ([6|), (O, and using the equation of state 



p/Pm, = and the conservation equation pm = poma 



(EOS) for matter 

equation in BD theory can also be expressed as [52 1 



-3(l+u),„) 



POmO. 



the Friedmann 



(10) 



with n = ^"^^^ ^' ^ constant parameters. Comparing Eqs. © and ([TU]). one has A = H'^Vlf^m, and 



the second term plays the role of potential. Thus Friedmann equation is written as, 

= i/o^f7o™a"(^+"' + Ba'\ 



(11) 



If interpret the constant parameter B in the second term as the current value of energy density for an " new" component 
including the BD potential, the above equation becomes 



^0 'Om" 



H^nQxa" 



(12) 
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with writing B = Hq^Iqx- According to Eq. p2|) and considering the index n = — 3(1 + Wj;), it is obtained that the 
equation of state is expressed as, Wx = —I — for this component. 




FIG. 1: The 2-D contours with la, 2a confidence levels and 1-D distribution of model parameters in the fiat BD (left) and CC 
(right) model. Solid lines are mean likelihoods of samples, and dotted lines are marginalized probabilities for ID distribution. 



In the following we apply the Markov Chain Monte Carlo (MCMC) method to investigate the global constraints on 
above BD scenario. For the used observational data, we consider 557 Union2 dataset of type supernovae la (SNIa) 
[ssj , 59 high-redshift Gamma- Ray Bursts (GRBs) data (mI, o bservational Hubble data (OHD) 55 1, X-ray gas mass 



fraction in cluster j56|, baryon acoustic oscillation (BAO) [57[, and cosmic microwave background (CMB) data [58j, 
and for the analysis method please see appendix. In our joint analysis, the MCMC code is based on the publicly 

6l|. The latter package is about the 



available CosmoMC package [59| and the modified CosmoMC package [56|, |60 
constraint code of X-ray cluster gas mass fraction. In the calculation the baryon matter density is taken to be varied 
with a tophat prior: Of,/i^ e [0.005,0.1]. In addition, we run 8 independent chains in the MCMC calculation, and to 
get the converged results we test the convergence of the chains by typically getting i? — 1 to be less than 0.03. The 
total is expressed as. 



Xtotal (Ps 

and the parameter vector reads 



XsNIa + XgRBs + XOHD + XcBF + XbAO + XcMBj 



Ps = {^Ibh^, ^ch?, a, n}. 



(13) 



(14) 



Here the expression of foi' each observation corresponds to Eqs.(jA5|, (jA6|) . (|A8j) . (jA12p . (jA20[) and (jA24|) . fif, and Vtc 
denote dimensionless energy density of baryon matter and dark matter, respectively. Based on the basic cosmological 
parameters Ps we can also obtain the derived parameters riom = + i^c, = 1 ^ ^^Om, and the Hubble constant 
Ha ~ lOO/i km-s ^^-Mpc^^. Using the currently observed data with the xlotai i'^ ^q. (fT3|). Fig. [T](left) plots the 2-D 
contours with la, 2a confidence levels and 1-D distribution of model parameters in the flat BD theory. Solid lines 
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are mean likelihoods of samples, and dotted lines are marginalized probabilities for ID distribution (for Gaussian 
distributions they should be the same) . And the calculation results are listed in table H] for the constraint on model 
parameters. In addition, Table HIl shows the values for the best-fit sample, and projections of the n-Dimcnsional la and 
2(7 confidence regions. The n-D limits give some idea of the range of the posterior, and are much more conservative 
than the marginalized limits [s^. Often, the best-fit results are recommended. According to the best fit values of 
model parameter, it is shown the value of equation of state is, = —1 — ^ ~ —1.034, which is in the phantom 
region with a small deviation. And according to the best fit values of a = 0.0046 and n = 0.102, the best fit value 
of coupling parameter, uj ~ -n(2+ay+2a(i-a) _ _/^Qir^ ii jg calculated, which is consistent with other results, such 
as w < —40 constrained from the growth rate data [62|, and u < —120 or lu > 97.8 constrained from the WMAP 
CMB and the SDSS BAO data 63]. In addition, comparing with the literature j^l, where the Brans-Dicke theory 
is constrained from the recent observational data with introducing the holographic dark energy in universe in order 
to obtain an accelerating universe, it can be found that our results in this paper have a larger value of parameter a. 
And for the constraint on the dimensionless matter density ilom in the BD theory, from Table HIl one can see that it is 
consistent with results from other analyses, such as th e const raints on several dynamical dark energy models and model 
independent scenarios by the recently observed data 
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67l |. where the best fit value of is about 0.27. At last, it 



may be said that the BD theory with a self-potential is slightly preferred by the current observational data because 
the quantity Xmin measures the goodness of model fit, i.e. the less value of Xmm' better model of agreeing with 
observations. Therefore, it seems that the BD theory with a self-potential is better than the cosmological constant to 
interpret an accelerating universe, and it is also significative to discuss the evolution of cosmological quantities in BD 
theory. 
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TABLE I: The means, standard deviations (the numerical results in brackets) and the marginalized limits for the model 
parameters from MCMC calculation, obtained by using SNIa Union2, GRBs, OHD, CBF, BAO, and CMB data. CC denotes 
the cosmological constant (CC) model in the Einstein's general relativity. 
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TABLE II: The maximum likelihood values Xmim Xmin / dof , the best fit model parameters, and the limits from the extremal 
values of N-dimensional distribution from the MCMC calculation by using SNIa Union2, GRBs, CBF, OHD, BAO, and CMB 
data, where dof is degree of freedom of model, and its value equals the number of observational data points minus the number 
of parameters. 



Considering the value of parameter a is not arbitrary, we calculate its values by using other experiment exploration 
on Newton s gravitational "constant" G. According to the independent observations of Hulse- Taylor binary pulsar 



B1913 -I- 16 



B^ls^, and astereoseismological data from the pulsating white dwarf star G117-B15A 



7C 



Q, it is 
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indicated by the current constraints on the variation of G, 



I g l< 10-"yr-i. 

For BD theory and the parameterized scenario (f) — (jjoo," , it corresponds to 

I ^ hi aH \< 10-"yr-i, 



(15) 



(16) 



I.e. 



I a |< X 10-"yr-i. 



(17) 



Considering the current value of Hubble constant h = 0.742lQ'Q3g [7^ and taking the best fit value of Hq = 74.2km-s 



'i-Mpc^i = 2.40 X lO^^^s"^ = 7.57 x 10""j/r"\ one obtains the bounds on a, 

I a |< 0.132124. 



(18) 



It is easy to see that the value of constraint on parameter a with 2(7 confidence level for BD theory is under this 
bound, i.e., it is lied in a physical significative region. 
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FIG. 2: The evolution of q(z) for BD and cosmological constant (CC) model. 
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TABLE III: The values of transition redshift zt, and current deceleration parameter qo from MCMC calculation, obtained by 
using SNIa Union2, GRBs, CBF, OHD, BAG, and CMB data. 



Furthermore, we investigate the evolution of deceleration parameter q{z) 



(l + z)-^-^ — 1 in this scenario. 



Considering the propagation of the errors for q{z) by the Fisher matrix analysis, the errors are evaluated by using the 
covariance matrix Cij of the fitting parameters (tsI \t\ , which is the inverse of the Fisher matrix and given by 



dpsidps^ 2 dps^dpsj 



(19) 
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where Ps is a set of parameters, and InL is the logarithmic hkelihood function. The errors on a function / = f{ps) in 
terms of the variables ps are calculated by [3, llal 

where m is the number of parameters, and / will be deceleration parameter q{z;ps-). The parameters respectively 
represent {flbh^ ,ilch^ ,a,n). In Fig. [5] (left) we plot the evolutions of q{z) with errors by 

Qiaiz) = q{z) \p^=p^ ±aq, (21) 

here ps are the best fit values of the constraint parameters. From this figure it can be seen that a universe from 
decelerated expansion to accelerated expansion is obtained. And the values of transition redshift zt and current 
decelerating parameter qo with la confidence level are, zt = 0.692^g'Q24 ^-^d qo = — 0.594^Q'Q3g, as shown in table IIIll 
For comparison we also consider the observed constraints on popular cosmological constant (CC), i.e. ACDM model 
in the Einstein's theory of GR (or dubbed as standard cosmology). The Friedmann equation for this case is expressed 
as 

— HQflomO, + iJo^oA- (22) 

The calculation results on cosmological parameters are listed in table |T] and |lT] for this model. In addition, for the 
evolution of deceleration parameter q{z) it is plotted in Fig. [2] (right), and the values of zt and qo are listed in table 

imi 



IV. jerk parameter geometrical diagnostic for Brans-Dicke theory 

In the following we use a new geometrical diagnostic method, jerk parameter j, to investigate discriminations 
between BD and CC model. The jerk parameter is defined by scale factor a and its third derivative 76-7^, 



The use of the cosmic jerk parameter provides more parameter space for geometrical studies, and transitions between 
phases of different cosmic acceleration are more naturally described by models incorporating a cosmic jerk. In addition, 
for a dark energy model Eq. ([23]) can be derived as 

9 3 ihde 

j ^ -I - -zWde^dei^ + Wde) + T^^de-rr^ (24) 
Z Z II 



where ^de denotes dimensionless energy density for dark energy. From Eq. (|24|) it is easy to see that, for fiat 
ACDM model {wde{z) ~ ^l)j it has a constant jerk with j{z) —1. Thus, it can provide us with a simple, 
convenient approach to distinguish and search departures for both cosmological dynamic and kinematical models 
from the cosmic concordance model, CC. Using Eq. (P^ . we plot the evolution of jerk parameter j{z) for BD model 
in Fig. [3] and compare it with cosmological constant. According to this figure we get the current values of jerk 
parameter, jo — — 1.002lg'Qg4. And it is shown that at la confidence level the BD theory can not be distinguished 
with cosmological constant model in standard cosmology according to the jerk parameter geometrical diagnostic. 
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FIG. 3: The evolution of jerk parameter j{z) for BD model. 
V. Conclusions 



Many observations and phenomena indicate that the gravitational theory of general relativity should be modified 



Scalar tensor (ST) theory as a modified theory of gravity has been widely studied 79|- 8^. Considering the ST theory 
can be well used to interpret the inflation problem, in this paper Brans-Dicke gravitational theory as an important 
one in ST theories is studied to interpret the late accelerating universe. In BD theory gravitational "constant" G is 
described by a function of scalar field which couples to the Ricci scalar R, so that the effective G generally varies in 
time. Concretely, with a parameterized scenario for BD scalar field 4> which plays the role of the variable gravitational 
constant G, we use the MCMC analysis method and the current observed data, including 557 Union2 SNIa, 59 
GRBs, OHD, cluster X-ray gas mass fraction, BAO and CMB data, to constrain the BD theory with a self-interacting 
potential. According to the constraint results on deceleration parameter q{z), the late accelerating universe can 
be obtained in this theory. And for the parameter a that is corresponded to the variable gravitational "constant" 
G in physics, its constraint value is consistent with the result of current experiment exploration, | a |< 0.132124. 
Considering the application of BD theory in cosmology, according to Eq. (fT2|) since the both energy densities are 
dynamical, not a constant as in cosmological constant model, it will not suffer the so-called cosmic coincidence 
problem. In addition, relative to other geometrical diagnostic methods, such as Om and statefinder parameters {r,s} 



which have been widely discussed in many models 89l - Qsl ] , in this paper we use a new geometrical diagnostic method — 
-jerk parameter to Brans-Dicke theory, and it is shown that it can not be distinguished with cosmological constant 
model. Obviously, the perturbation theory could distinguish this modified gravity theory of BD with the CC dark 
energy scenario in Einstein's theory of general relativity, which deserves to be studied in the future. 
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Appendix A: Observational data and cosmological constraint methods 

In this part we introduce the cosmological constraints methods and the current observed data used in this paper. 



1. Type la supernovae 

j 

We use the 557 SNIa Union2 dataset, which includes 557 SNIa 53]. Following 
sponding constraints by fitting the distance modulus as 

15 G 

^j,th{z) = 51ogio[DL(2)] + — logio — +^J.o, 

4 Go 



96h99||. one can obtain the corre- 



(Al) 



where Gq is the current value of Newton's constant G. In this expression Dl{z) is the Hubble- free luminosity distance 
HodL^z) / c, with Hq the Hubble constant, defined through the re-normalized quantity h as Hq = lOO/i km s^^Mpc^^, 
and 

, ■ C(l + Z) . r— dz' 

dL[z) = — j==^smn[y'\Uk\ 



Hi^T 



MO = BlogioC-r^) + 25 = 42.38 - 5logioh. 
Mpc 



where s^nnn{^J\Vlk\x) respectively denotes sin(-\/|f2fc|x), -^/jrifelx, sinh(-\/|f2fe|x) for fifc < 0, 51^ = and J7fe > 0. 
Additionally, the observed distance moduli ^iohs{zi) of SNIa at Zi is 



(A2) 



where M is their absolute magnitudes. 

For using SNIa data, theoretical model parameters Ps can be determined by a likelihood analysis, based on the 
calculation of 



X 



itJ-obsiZi) - P-thiPs, Zi)} 



E 



{5\0gj^Q[DL{Ps,Zi)] - TTlobsiZt) + M'Y 



(A3) 



SNIa * SNIa 

where M' = /io + M is a nuisance parameter which include s the abs olute magnitude and the parameter h. The 
nuisance parameter M' can be marginalized over analytically 



IOC -1061 as 



/ + 00 
exp 
-OO 



-\x\Ps.M') 



dM\ 



resulting to 



with 



2 . . (C 
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{51ogio[-DL(Ps,^i)] -"^ofcs(^^)} 



SNIa 



51ogio[-DL(Ps,2:i)] -mobsizi) 



SNIa 



SNIa 
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Relation (|A3[) has a minimum at the nuisance parameter value M' = B /C , which contains information of the values 
of h and M . Therefore, one can extract the values of h and M provided the knowledge of one of them. Finally, note 
that the expression 

xlNIaiPs)=A-{ByCl 



96, 



97 



IOC -1031, and thus in this 



which coincides to (|A4p up to a constant, is often used in the likelihood analysis 
case the results will not be affected by a flat M' distribution. For minimizing X%Nia^Ps) to perform a constraint on 
cosmological parameters, it is equivalent to maximizing the likelihood 



L{ps) oc exp[ 



(A5) 



2. high-redshift Gamma-Ray Bursts data 

The GRBs data can be observed at higher redshift than SNIa. The currently observed reshift range of GRBs is at 
0-^ ^ z < 9. Therefore, the GRBs data can be viewed as an excellent complement to SNIa data and would provide 
more information at high redshift. When several empirical relations of the GRBs are proposed, these indicators have 
motivated the authors make use of the GRBs as cosmological standard candles at high redshift. However, the fact that 
there are not sufficient low reshift GRBs leads that the calibration of GRB relations is dependent on the cosmological 
model, namely, the circularity problem. One of methods to solve the circularity problem is the calibrat ion o f GRB 
relations are performed by the use of a sample of SNIa at low redshift in the cosmology- independent way[l07|. Here, 
the GRBs data we used consists of 59 GRB samples with a redshift range of 1.4 < z < 9 obtained in [54|. These 59 
GRBs are calibrated by utilizing the new ly re leased 557 Uion2 SNIa and the isotropic energy-peak spectral energy 
{Eiso- Ep,i) relation (i.e. Amati relation) [108|. 

The Xgrbs takes the same form as x'sNia 



59 



[^J■obs{zi - Hth{zi;ps,lJ.o)f 



(A6) 



The same method are used to deal with the nuisance parameter /^o a-s shown in the description of x'sNia above. 



Observational Hubble data 
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The observational Hubble data [l09| are based on differential ages of the galaxies . In [llo| , Jimenez et al. obtained 



and used it to constrain the 



an independent estimate for the Hubble parameter using the method developed in 
cosmological models. The Hubble parameter depending on the differential ages as a function of redshift z can be 
written in the form of 



H{z) = 



1 dz 
1 + z dt' 



(A7) 



So, once dz / dt is known, H{z) is obt ained directly. By using the differ ential ages of passively-evolving galaxies from the 
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llSf . Simon et al. obtained several values oi H(z) at 



Gemini Deep Deep Survey (GDDS) [ll2!| and archival data 
different redshift 55j. The twelve observational Hubble data (redshift interval < z < 1.8) from 
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l3 



are list 
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0.1 


0.17 


0.27 


0.4 


0.48 


0.88 


0.9 


1.30 


1.43 


1.53 


1.75 


H{z) (km s-^Mpc-i) 


74.2 


69 


83 


77 


95 


97 


90 


117 


168 


177 


140 


202 


la uncertainty 


±3.6 


±12 


±8 


±14 


±17 


±60 


±40 


±23 


±17 


±18 


±14 


±40 



TABLE IV: The observational H(z) data |72l. Ill9l. ll2C | 



in Table iTVl In addition, in [121" the authors take the BAO scale as a standard ruler in the radial direction, and obtain 
three more additional data: H{z = 0.24) 79.69 ± 2.32, H{z = 0.34) = 83.8 ± 2.96, and H{z = 0.43) = 86. 45 ±3.27 
The best fit values of the model parameters from observational Hubble data are determined by minimizing 
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Xohd{Ho,Ps 



[Hth{Ho,Ps]Zi) - ffobs(Zi)]' 



<y^z,) 



(A8) 



where Htu is the predicted value for the Hubble parameter, Hobs is the observed value, cr(zi) is the standard deviation 
measurement uncertainty, and the summation is over the 15 observational Hubble data points at redshifts Zi. 



4. The X-ray gas mass fraction 

The X-ray gas mass fraction, fgas, is defined as the ratio of the X-ray gas mass to the total mass of a cluster, which 
is approximately independent on the redshift for the hot (fcT > bkeV)^ dynamically relaxed clusters at the radii larger 
than the innermost core r25oo- As inspected in [56|, the ACDM model is very favored and has been chosen as the 



reference cosmology. The model fitted to the reference ACDM data is presented as [5 



KA-fb{z) ( n 



1 + s{z) \il,o 



Da{z) 



(A9) 



where where D'^^'^'^ [z) and Da{z) denote respectively the proper angular diameter distance in the ACDM reference 
cosmology and the current constraint model. A is the angular correction factor, which is caused by the change in 
angle for the current test model 6*2500 in comparison with that of the reference cosmology 6*25^0*^^ • 



A = 



I 02500 ) ^ \[H{z)Da{z)]^^^'' ) ' ^ ^ 



here, the index 77 is the slope of the f gas{i' / f2boo) data within the radius r25oo, with the best-fit average value 
77 = 0.214 ± 0.022 56^. And the proper (not comoving) angular diameter distance is given by 



Da{z)^- ^^sinn[/|f7^ / — -]. (All) 

{l + z)y/\yi^\ Jo H{z') 



It is clear that this quantity is related with d^i^z) by 

dL{z) 



Da{z) 



{l + zf 

In equation (|A9|) . the parameter 7 denotes permissible departures from the assumption of hydrostatic equilibrium, 
due to non-thermal pressure support; the bias factor b{z) — feo(l ± otbz) accounts for uncertainties in the cluster 
depletion factor; s{z) — so(l ± Ugz) accounts for uncertainties of the baryonic mass fraction in stars and a Gaussian 
prior for so is employed, with so = (0.16 ± Q.Qh)h^Q the factor K is used to describe the combined effects of the 
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residual uncertainties, such as the instrumental calibration and certain X-ray modelling issues, and a Gaussian prior 
for the 'calibration' factor is con sidered b y K = 1.0 ± 0.1 ^ 
Following the method in Ref. 



56 



125 



and adopting the updated 42 observational fgas data in Ref. 



561, the best 



fit values of the model parameters for the X-ray gas mass fraction analysis are determined by minimizing, 



XCBF 



N 



[f^.?s''^'i^^)-fsas{Zi)]^ , (sq "0.16)2 , (X-l.O)^ - 0.214)^ 



0.00162 



0.012 



0.0222 



(A12) 



where crfg^^{zi) is the statistical uncertainties (Table 3 of 56]). As pointed out in :5(|, the acquiescent systematic 
uncertainties have been considered according to the parameters i.e. 77, b{z), s{z) and K. 



5. Baryon acoustic oscillation 



The baryon acoustic oscillations are detected in the clustering of the combined 2dFGRS and SDSS main galaxy 
samples, which measure the distance-redshift relation at zbao = 0.2 and zbao = 0.35. The observed scale of the 
BAO calculated from these s ampl es, are analyzed using estimates of the correlated errors to constrain the form of the 



distance measure Dv{z) [571 Il26 



Dy{z)^[{l + zyD'Jz) 



dz 



2li 



(A13) 



' H{z;ps)' E{z;ps) Jq E{z';ps)' 

In this expression E(z;ps) — H{z;ps)/Ho, Da{z) is the proper (not comoving) angular diameter distance, which has 
the following relation with dL(z) 

DJz) - ^^^^^ 

The peak positions of the BAO depend on the ratio of Dv{z) to the sound horizon size at th e drag epoch (where 
baryons were released from photons) z^, which can be obtained by using a fitting formula 129|: 



Zd 



1 + 0.659(fiom/l2)0-828 



[l + 6i(f]6/l2 



with 



61 = 0.313(r!om/'-')"°-^'^[l + 0.607(17o™/i2)0-674] 

62 = 0.238(r!o™/(,2)0-223^ 



(A14) 

(A15) 
(A16) 



In this paper, we use the data of rs{zd)/ Dv{z) extracted from the Sloan Digitial Sky Survey (SDSS) and the Two 
Degree Field Galaxy Redshift Survey (2dFGRS) 126l4l28l |. which are listed in Table IVl where rs{z) is the comoving 
sound horizon size 



rs[z) =c / 



csda r- ^ c 
c I -TT— — CI dz 



H{z) 



c 

TsJo 



1/(1 + 2) 



da 



where Cs is the sound speed of the photon— baryon fluid 



130 



13l|: 



c„ - = 3 



4 ^ pb{z) 
3 p^{z) 



(A17) 



(A18) 
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z 


r,(2d)/Dv(2) 


0.2 


0.1905 ±0.0061 


0.35 


0.1097 ±0.0036 



TABLE V: The observational rs{zi)IDv{z) data 



E3^ 



30124.1 -17226.9 
-17226.9 86976.6 



(A19) 



and here VL^ = 2.469 x lO^^/i"^ for Tcmb^ 2.12bK. 

Using the data of BAO in Table IVl and the inverse covariance matrix in [57[: 

the XBAoiPs) is given as 

xlAoiPs)=X'V-'X, (A20) 
where X is a column vector formed from the values of theory minus the corresponding observational data, with 



X 



Dv(0.2) 
Dv(0.35) 



0.1905 
-0.1097 



(A21) 



and X* denotes its transpose. 



6. Cosmic microwave background 

The CMB shift parameter R is provided by 132 1 



Hf)dz 



R = Jno„,H§{l + z,)Da{z,)/c = / r, (A22) 

Jo ^ 'Ps) 

here, the redshift (the decoupling epoch of photons) is obtained using the fitting function |l33t 

= 1048 [1 + 0.00124(r2f,/i2)-"-738] [l + gi{n„,h^)3^] , 

where the functions gi and g2 read 

gi = 0.0783(f}6/i')-°'''(l±39.5(r!b/i2)0.763)-i^ 

52 = 0.560 (1 ± 21. 1(r2b/i2)i si)"\ 

In addition, the acoustic scale is related to a distance ratio, Dji{z)/rg{z), and at decoupling epoch it is defined as 

Z^^(l + z,)^^^fM, (A23) 

where Eq. (jA23P arises a factor 1 + z*, because Da(z) is the proper (physical) angular diameter distance, whereas rs(z*) 
is the comoving sound horizon. Using the data of Ia,R,z^ in [58] and their covariance matrix of [^^(2*), -R(2:*), z*] 
(please see table IVII and IVII|) . we can calculate the likelihood L as Xcmb = — 21nL: 

xImb = Ad,[Cov'\d,,dj)[Ad,]% (A24) 



14 





7- year maximum likelihood 


error, a 


Ia{z,) 


302.09 


0.76 


R{z,) 


1.725 


0.018 


z* 


1091.3 


0.91 



TABLE VI: The values of Ia{z*), R{z*), and z, from 7-year WMAP data. 





Ia{z*) 


R(z.) 


2* 


lAiz,) 


2.305 


29.698 


-1.333 






6825.270 


-113.180 








3.414 



TABLE VII: The inverse covariance matrix of Ia{z*), R{zt), and z* from 7-year WMAP data, 
where Adi = di — df°'^°' is a row vector, and di = [Ia, R,z^). 
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